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Abstract 

We study the dynamics of inertial particles in two dimensional incompressible flows. The Maxey- 
Riley equation describing the motion of inertial particles is used to construct a four dimensional 
dissipative bailout embedding map. This map models the dynamics of the inertial particles while 
the base flow is represented by a 2 — d area preserving map. The dynamics of particles heavier than 
the fluid, the aerosols, as well as that of bubbles, particles lighter than the fluid, can be classified 
into 3 main dynamical regimes - periodic orbits, chaotic structures and mixed regions. A phase 
diagram in the parameter space is constructed with the Lyapunov characteristic exponents of the 
4 — d map in which these dynamical regimes are distinctly identified. The embedding map can 
target periodic orbits, as well as chaotic structures, in both the aerosol and bubble regimes, at 
suitable values of the dissipation parameter. 
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I. INTRODUCTION 



The motion of inertial particles in fluids is governed by dynamical equations which display 
rich and complex behaviour. The inertial particles have density that is different from that 
of the fluid in which they are immersed, in contrast to the neutral particles which have 
the same density of the fluid. The dynamics of the inertial particles deserve attention from 
both the point of view of fundamental physics as they exhibit complex behaviour, as well 

situations e.g. the transport of pollutants in 



as from that of their applicability to practica 
the atmosphere and plankton in oceans 1, [2I, 



If impurities in fluids are modeled by small spherical tracers, the Lagrangian dynamics of 
such small spherical tracers in two dimensional incompressible fluid flows is described by the 
Maxey-Riley equations. These are further simplified under various approximations to give 
a set of minimal equations called the embedding equations where the fluid flow dynamics 
is embedded in a larger set of equations which include the differences between the particle 
and fluid velocities [l, 4]. Although the Lagrangian dynamics of the underlying fluid flow is 
incompressible, the particle motion is compressible 7|], and has regions of contraction and 
expansion. The density grows in the former giving rise to clusters and falls in the latter 
giving rise to voids. The properties of the base flow have important consequences for the 
transport and mixing of particles. Map analogs of the embedding equations have also been 
constructed for cases where the fluid dynamics is modeled by area-preserving maps which 
essentially retain the qualitative features of the flow @, 9|. The embedded dynamics in both 
cases is dissipative in nature. 

Further complexity is added to the dynamics by the density difference between the parti- 
cles and the fluid. In the case of two-dimensional chaotic flows, it has been observed earlier 
that particles with density higher than the base flow, the aerosols, tend to migrate away 
from the KAM islands, while particles lighter than the fluid, the bubbles, display the op- 
posite tendency pi]]. Neutrally buoyant particles also showed a similar result, wherein the 
particles settled into the KAM islands. Our study indicates that in addition to the density 
difference, the dissipation parameter of the system also has an important role to play in the 
dynamic behaviour of the aerosols and the bubbles. 

In Sec. II of this paper, the Maxey-Riley equation describing the motion of inertial 
particles is simplified to get a minimal equation of motion. A map analog of this minimal 
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equation called the embedding map, is obtained to model the particle dynamics in discrete 
time. The fluid dynamics is represented by a 2 — base map, which we choose to be the 
standard map. The embedding map is four dimensional and dissipative, while the base map 
is two dimensional and area-preserving. Unlike the earlier results mentioned above, viz. that 
the bubbles tended to form structures in the KAM islands[10] and the aerosols were pushed 
away, we found that structures form for both bubbles and aerosols in certain parameter 
regimes due to the role of the dissipation parameter. Both the aerosol and bubble regimes 
in the phase diagram show regions where periodic orbits as well as chaotic structures can 
be seen in the phase space plots. In Sec. Ill we obtain the phase diagram of the system in 
the a — 7 space where a is the mass ratio parameter, and 7 is the dissipation parameter. 
The phase diagram shows rich structure in the a < 1 aerosol regime, as well as the a > 1 
bubble regime. In addition to these, fully or partially mixed regimes can also be seen in 
both the aerosol and bubble cases. Thus the dynamic behaviour of the inertial particles 
can be of three major types. The Lyapunov exponents of the four dimensional map is 
used as the characterizers in demarcating these three major regimes. Our results can have 
implications for practical problems such as the dispersion of pollutants by atmospheric flows, 
and catalytic chemical reactions. 



II. THE EMBEDDING EQUATION 



The motion of small spherical particles in an incompressible fluid has been studied ex- 
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Here v represents the particle velocity, u the fluid velocity, pp the density of the particle, 
Pf the density of the fluid, and z/, a, g represent the kinematic viscosity of the fluid, the 
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radius of the particle and the acceleration due to gravity respectively. 

The first term in the right of Eqn. [1] represents the force exerted by the undisturbed 
flow on the particle, the second term represents the buoyancy, the third term represents the 
Stokes drag, the fourth term represents the added mass, and the last the Basset-Bossinesq 
history force term. The derivative Du/Dt is taken along the path of the fluid element, 
Du/Dt = du/dt + (u ■ V)u. The derivative du/dt, is taken along the trajectory of the 
particle du/dt = du/dt + {v ■ V)u. 



In deriving Eqn. [T] it was assumed [3|, [llj that the particle radius, and the Reynolds 



number are small as well as the velocity gradients around the particle. It is also assumed 
that the initial velocities of the particle and the fluid are same. A full review of the problem 



can be found in Ref. 



12| 



Now Eqn. [T] can be simplified through the following arguments. The Faxen correction 
a^V^u is of the magnitude 0{o?u/ L), and from the assumption, a « L, this term's con- 
tribution becomes negligible and can be excluded from the equation. The Basset history 
force term which takes into account viscous memory effects becomes less significant and can 
be dropped, as the particle size is sufficiently small and the concentration of particles is 
sufficiently low, that they do not modify the fiow field or interact with each other 
Under the low Reynolds numbe r ap proximation , both the derivatives Du/Dt and dv/dt 



will approximately be the same 



ll|. We assume the buoyancy effects to be negligible. We 



thus arrive at the following equation. 



dv du '^vpf , Pf fdv du\ . . 

This can be easily written into a non-dimensional form by proper rescaling of length, 
time and velocity by scale factors L,T = L/U ,and U respectively, to obtain. 



dv du 2 f 9a \ . . . . 

where a = 3p//(2pp + P/) is the mass ratio parameter. Re = {UL)/u, is the Reynolds 
number of the fluid. By using the particle Stokes number St = {2/9)a?Re and deflning the 
dissipation parameter as 7 = (2a/3S't), the flnal simplifled equation takes the form. 
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FIG. 1: The phase space plot of the standard map for K = 2.0. Random initial conditions 
distributed uniformly in the phase space were chosen. 



dv dvL , . 

The map analog of Eqn. H] can be constructed. Let the base flow be represented by an 
area preserving map M with the map evolution equation x„+i = M(x„). As the particle 
dynamics in the fluid flow is compressible, it can be represented by a dissipative map. The 
map analog of Eqn. IHhas been chosen in |l| to be, 



x„+2 - M(x„+i) = e ^{ayin+i - M(x„)). (5) 

This equation can be rewritten by introducing a new variable 5, which defines the de- 
tachment of the particle from the local fluid parcel as, 

X„+i = M(x„) + 5n 

5n+i = e~^[axn+i - M(x„)]. (6) 

This is the bailout embedding map. When the detachment measured by the 6 is nonzero, 
the particle is said to have bailed out of the fluid trajectory. In the limit 7 ^ oo and 
a = 1, (5 — i> and the fluid dynamics is recovered. So the fluid dynamics is embedded in the 



particle's dynamics and is recovered under appropriate limit. This map is dissipative with 
the phase space contraction rate being e~^". 

We choose the standard map 13(], to be our base map M as it is a prototypical area- 
preserving system and is widely studied in a variety of problems of both theoretical and 
experimental interest. The standard map serves as a test bed for various models regarding 
transport 1^. The map equation is given by, 



x„+i = x„ + y„+i {Mod 1) 

K 

Yn+i = y„ + — sin(27rx„) (Modi). (7) 
Ztc 

The parameter K controls the chaoticity of the map. This map is taken to be the base 
map M(x^) in Eqn. (3). The phase space of the base standard map for K = 2.0 is plotted 
in Fig. [T]in the region x G [—1, l],y G [—1, 1]. The particle dynamics is governed by the 4 
dimensional map , 

K 

Xn+i = Xn + — sin(27ry„) + 6^ 
Zn 

K 

Un+i = Xn + yn + — sin(27ry„) + SI 

= e-"'[aXn+l - {Xn+1 - 6^)] 

Sl+i = e-^[ayn+i - {Vn+i - S^)] (8) 

It is clear that this 4-dimensional map is invertible and dissipative. The embedding map 
has 3 parameters K, a, and 7. We plot the phase space portrait of the system evolved with 
random initial conditions at i^' = 2.0. The phase space plot of the standard map at the 
same value of K in Fig. [T] shows islands and chaotic regions separated by invariant tori 
which constitute barriers to transport. We plot the embedding map for the aerosol and the 
bubble cases, in Fig. [2l Unlike earlier studies which claim that only the aerosols are able to 
breach the invariant tori, we observe that that both the aerosols and bubbles have broken 
the barrier posed by the invariant curve and have targeted the islands forming periodic 
trajectories and structures. Thus, it is clear that the dissipation parameter 7 plays a vital 
role in the behaviour of both kinds of particles, the aerosols and the bubbles. Depending on 
particular values of 7 and a both these particles can form clusters or undergo mixing in the 
phase space. Therefore in order to identify the regimes in which clustering or mixing can 
take place, it is necessary to study the full phase diagram in the parameter space. 
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FIG. 2: Phase space plots showing four different dynamical regimes of the embedding map. Here 
K = 2.0 is fixed, (a) periodic trajectories (a = 1.5,7 = 0.7), (b) chaotic structure (a = 0.2,7 = 
0.9), (c) full mixing (a = 2.0,7 = 0.9)and (d) chaotic structure in a mixing background (a = 
0.7,7 = 0.18), 



III. THE PHASE DIAGRAM 



The phase diagram of the system can be constructed using the following criteria. The 
regimes where the largest Lyapunov characteristic exponent (LCE) Xmax > are regions of 
chaos and the regimes where Xmax < are regions of periodic behaviour. This property is 
used to identify the periodic regimes (marked in blue) in the phase diagram (see Fig. [3l). A 
typical phase space plot in the periodic regime is shown in Fig. [2] (a). For the calculation 
of the Lyapunov exponent, the first 5000 iterates are taken as transients and the next 1000 
iterates are stored, for 100 random initial conditions, uniformly spread in the phase space. 



7 



We use the Gram-Schmidt reorthonormahzation (GSR) procedure for the calculation of all 
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the Lyapunov exponents of the embedded map 

In the chaotic regimes, further classification can be carried out based on a scheme which 
uses the signs of the other Lyapunov exponents of the system and the occupation densities of 
the trajectories in the phase space. Fig. H] shows all the Lyapunov exponents of the system 
for the dissipation 7 = 0.7 and a = — 3. Chaotic structures in the phase space are found 
in the regimes where A2,3,4 < while the largest LCE Ai (i.e.,Amax) is greater than 0. These 
are marked in green in the aerosol regions and yellow in the bubble regions in the phase 
diagram. For the regimes where both Ai and A2 greater than 0, mixing is seen in the phase 
space. 

Fig. [2] shows the phase space plots seen in three different dynamical regimes. Fig. 
[2t^a) shows the periodic orbits in the bubble regimes [a = 1.5,7 = 0.7). Here, both Ai 
and A2 are less than 0. Chaotic structures are shown in Fig. [2](b) in the aerosol regime 
{a = 0.2,7 = 0.9), where Ai > and A2 < 0. The Figs. [2] (c),(d) show the mixing regimes 
for the bubble (a = 2,7 = 0.9) and the aerosol (a = 0.7,7 = 0.18) regimes respectively. In 
both the cases for the mixing regime, Ai and A2 are greater than 0. We also note that the 
Lyapunov exponents A3 and A4 are less than in all the three regimes. 

We see that the criterion based on the signs of Lyapunov exponents cannot distinguish 
between the phase space plots showing mixing regime (Fig. [2] (c)) and the chaotic structure 
in a mixing background (Fig. [2](d)). Therefore we use the occupation densities to distinguish 
between these two cases. 

The occupation densities in phase space can be calculated by covering the phase space 
by a 100 x 100 mesh, and counting the number of boxes accessed by the iterates. The initial 
conditions and the transients are same as those that were used for calculating the Lyapunov 
exponent. The regions of parameter space where the iterates access more than 99% of the 
phase space grid turn out to be fully mixing regions (marked in red in Fig. [3]), with no 
remnant of any chaotic structure seen in the phase space plot, and have been marked with 
an M. The phase space plot of a mixing regime is seen in Fig. O The regions where the 
iterates access more than 90% and less than 99% of the boxes in the phase space are seen 
to contain a chaotic structure in a mixing background (marked with a '*' - light blue in Fig. 
[3]) . Fig. [2] (d) shows the phase space plot of a chaotic structure in a mixing background. 

Periodic structures are seen inside tongues in the parameter space in the aerosol as well 
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FIG. 3: (Color Online) The phase diagram. The regimes with periodic structures are marked with 
P (blue). The regimes with chaotic structures are marked with C (green for aerosols and yellow for 
bubbles). The regimes with structures in a mixing background are marked with a '*'(light blue). 
The regimes with full mixing are marked with M (red). 

as bubble regimes. Thus, the clustering or preferential concentration of inertial particles can 
be seen in both the aerosol and bubble regimes. Chaotic structures can be seen outside the 
tongues in the aerosol regime and in the border regions of the tongues in the bubble regime. 
The aerosol part (0 < a < 1) of the phase diagram has a small region of mixing while the 
)ubble part (1 < a < 3) is dominated by the mixing regime. In recent work we reported 
16(1 that crisis induced intermittency is seen at some parameters in the aerosol regime in 
the phase diagram. Here the characteristic time between bursts scales algebraically with a 
function of the dissipation parameter 7, with the power law exponent to he (3 = —1/3. 



IV. CONCLUSION 



The present work discusses the Lagrangian dynamics of passive scalars of finite size in 
an incompressible two-dimensional flow for situations where the particle density can differ 
from that of the fluid. The motion of the advected particles is represented by a dissipative 
embedding map with the area preserving standard map as the base map. 
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FIG. 4: (Color Online) Four Lyapunov exponents of the embedding standard map; K = 2.0 and 
7 = 0.7. Intervals for which Ai^2,3,4 < have periodic orbits, and X2,3a < have chaotic structures 
and Ai^2 > have mixing regimes in the phase space. 

The embedded standard map has three parameters, the nonlinearity K, the dissipation 
parameter 7 and the mass ratio parameter a which measures the extent to which the density 
of the particles differs from that of the fluid. The phase diagram of the system in the 
a — 7 parameter space shows very rich structure with three types of dynamical behaviours - 
periodic orbits, chaotic structures, and mixed regimes which can be partially or fully mixed. 
The Lyapunov characteristic exponents (LCE) of the A — d map are used to demarcate these 
dynamical regimes in the phase diagram. We observed that both aerols and bubbles can 
breach the invariant curves of the base map and form clusters, depending on the dissipation 
parameter of the system. 

Our study can have implications for the preferential concentration of inertial particles in 
flows, as well as for the targeting of periodic structures. Earlier results for such systems 
indicated that bubbles could breach elliptic islands and target structures, whereas aerosols 
could not. Our results indicate that both aerosols and bubbles can breach the invariant 
curves of the base map and form clusters, depending on the dissipation parameter of the 
system. Mixing regimes can also exist for both aerosols and bubbles. Thus the dissipation 
parameter 7 plays as crucial role in clustering and mixing as does the density differential 
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parameter a, and the examination of the full phase diagram is necessary to draw conclusions 
about the parameter regimes where targeting and breaching can occur. These results can 
have implications for the dynamics of impurities in diverse application contexts, e.g. that 
of the dispersion of pollutants in the atmosphere, flows with suspended micro-structures, 
coagulation of material particles in flows, catalytic chemical reactions and the plankton 
population in oceans. 
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